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Abstract. We present a study of the two species totally asymmetric diffusion model using 
the Bethe ansatz. The Hamiltonian has U q (SU(3)) symmetry. We derive the nested Bethe 
ansatz equations and obtain the dynamical critical exponent from the finite-size scaling 
properties of the eigenvalue with the smallest real part. The dynamical critical exponent 
is | which is the exponent corresponding to KPZ growth in the single species asymmetric 
diffusion model. 
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1 Introduction 

The single species asymmetric diffusion process has attracted and continues to attract a lot of 
interest. It is a driven diffusive system describing a stochastic movement of hard-core particles 
in one dimension, where movements to the left and right take place with different probabilities. 
There are many interesting applications to shock formation [H [21 [3] , traffic flow [3] , biopoly- 
mers [5j EJ [7] and other driven diffusive systems^). Since the underlying quantum spin chain is 
integrable, the powerful method of Bethe ansatz is applicable and was widely used in the past 
to obtain analytical results (see e.g. [8j EJ [TOJ [TT| I12j). The two-species asymmetric diffusion 
problem however has been less investigated. While the stationary state has been studied in some 
detail [I3j [TU HU CEU [17], much less is known about the full dynamics of the system. 

Some results for the dynamical phase diagram were obtained from the matrix product 
ansatz [18] (however these are restricted to certain regions of the full parameter space) and 
from numerical studies of the model [19] . 

Recently, the multi-species asymmetric diffusion model with open boundaries has sparked 
interest [20j (2TJ [22] . This case can also be tackled using the matrix product ansatz. 

In [23] the dynamics of a particular choice of the model is studied where two species hop 
in opposite directions on a ring, the diffusion constant being different from the passing rate 
of particles of different species |24j . This model corresponds directly to two coupled Burgers 
equations or to a Burgers equation coupled to a diffusion equation. The dynamical scaling 
properties are calculated from the time-evolution of the two-point correlation functions, based 



*This paper is a contribution to the Proceedings of the XVIIIth International Colloquium on Integrable Sys- 
tems and Quantum Symmetries (June 18-20, 2009, Prague, Czech Republic). The full collection is available at 
|http://www.emis.de/journals/SIGMA/IS QS2009.html 

i We are not able to give an exhaustive list of the vast literature about this subject here. 
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on a numerical study. The dynamical critical exponent is consistent with a subtle double Kardar- 
Parisi-Zhang type scaling (in the sense of factorization) . 

An analytical study of the dynamical scaling of the multi-species asymmetric exclusion pro- 
cess was described by [25] who found a dynamical critical exponent of | from a Bethe ansatz 
calculation if the diffusion was asymmetric and a critical exponent of 2 for symmetric diffusion. 

In the present paper, we investigate the dynamical critical exponent for the totally asymmetric 
diffusion process with equal rates which is a related but dynamically different model that will 
be precisely defined in Section [2j We will derive the Bethe ansatz equations starting from the R- 
matrix of the totally asymmetric diffusion model, calculate the energy gap and from a finite-size 
scaling analysis extract the critical exponent. 

1.1 Discussion of the literature about the two-species 
asymmetric diffusion model 

We will briefly survey the literature on the two-species asymmetric diffusion model. Although 
there seem to be several approaches building on previous work [2SJ [2UJ [27J ES] by trying to 
extract an appropriate limit to totally asymmetric diffusion, the details of each such strategy 
are not obvious, however, and one runs into serious difficulties and subtleties as we will explain. 
We found that instead of trying to tweak the calculations in an ad-hoc manner, it is clearer 
and simpler to start from scratch and derive the Bethe ansatz equations directly for the totally 
asymmetric case. With hindsight, our calculation provides a firm ground for future comparison 
which might also facilitate to formulate the correct limiting procedure. 

The Bethe ansatz equations for the two-species asymmetric diffusion model were put forth 
in an article by Dahmen j26j . They depend on a parameter 7 that is related to the asymmetric 

hopping rates by exp(7) = q = w p^-, where Tl and Tr describe the hopping rate to the left 

y R 

and right, respectively. It is not at all clear how to take the limit of one of these hopping rates 
going to zero directly in the Bethe ansatz equations. Since in the limit q — > all the ratios tend 
to one, and even the usual rescaling of the roots by a factor of q does not work - the resulting 
equations either diverge or become trivial. 

In the paper by Popkov at al. j27j Bethe ansatz equations for three-species models with 
different hopping rates for all six possible particles exchanges are derived from the matrix 
product ansatz. One of the cases they discuss is a particular asymmetric exclusion process 
where the hopping rates are coupled as gBo = 9A0 + 9AB where gso, gAO and gAB correspond to 
the processes A + ^- + A, -B + — > + B and A + B — >• B + A, respectively (and the three 
other rates are set to zero). In this paper we will choose all three non-zero rates to be the same, 
so the model considered here is different. 

Although the way Cantini [28] derives the Bethe ansatz equations for the case of asymmetric 
diffusion seems to be the most amenable for taking the limit to our model, however, he has 
introduced phase factors e^ 10 , e^ 20 , e Vl2 that remain in the Bethe ansatz equations, and it is not 
clear which value they should take for the case of totally asymmetric diffusion. 

In the article by Arita et al. |25] the Bethe ansatz equations for the multi-species asymmetric 
diffusion model were derived. However, the Bethe ansatz equations for the case of totally 
asymmetric diffusion are impossible to obtain from this calculation in a straightforward manner. 
Although it is possible to take the limit where one of the hopping rates is set equal to zero at 
the starting point, the i?-matrix, in the subsequent calculations the authors often multiply by or 
divide by the parameter that would become zero in the totally asymmetric limit. So taking the 
limit in the final equation does not lead to reasonable equations, and would have to do a step by 
step analysis of the derivation and make the necessary alterations at each step to avoid division 
by or multiplication by zero. At this point, we found it simpler and less ambiguous to just derive 
the equations for our case directly. 
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Figure 1. Example of totally asymmetric diffusion on a ring. Black dots correspond to A particles, grey 
dots to B particles and open dots to vacancies. 

2 Two-species totally asymmetric diffusion 

2.1 Asymmetric diffusion 

The two-species asymmetric diffusion model consists of two species of particles, A and B, 
diffusing asymmetrically in one dimension. The following processes take place: 

A + -)■ + A with rate F R , + A -)• A + with rate T L , 

B + 0-^0 + 5 with rate T R , + B ->• B + with rate T L , 

A + B ->• B + A with rate F R , B + A-> A + B with rate T L . 

2.2 Totally asymmetric diffusion 

In the totally asymmetric case, Tl = 0, so particles A do not diffuse to the left since the 
interchange B + A — > A + B is blocked. This leads to a different dynamics for the model that 
will be analyzed in this article. 

Assuming periodic boundary conditions, i.e. on the ring Z/LZ (where L is the number of 
discrete sites of that ring), the model is shown qualitatively in Fig. [TJ Diffusion to the right 
corresponds to clockwise motion around the circle, diffusion to the left to counterclockwise 
motion. Black dots represent particles A, grey dots particles B and open dots vacancies. The 
arrows indicate processes that are still allowed of Yl = 0, the blocked arrow indicates the process 
B + A — > A + B that is forbidden in the totally asymmetric case. 

Particles A do not see any difference between vacant sites and particles B, whereas particles B 
trying to diffuse to the right are blocked by particles A. Therefore, we call particles A "first- 
class", and particles B "second-class" particles. There is an interesting connection between 
second-class particles and the study of current fluctuations [29] . 

2.3 Master equation and Hamiltonian 

The dynamics of the asymmetric diffusion model is described by a master equation for the pro- 
bability distribution p t (77) of the lattice configuration rj(t) at time t. If we denote a configuration 
by rj and the jump rates T R and Tl by c(j,j + l,rf), when interchanging the particles and/or 
vacancies of the configuration 77 on sites j and j + 1, the master equation reads: 

j+Ptiv) = J2 K 7 '-? + 1 ^ 7 i n+1 )pt( 1 i n+1 ) - c(i, i + 1,77)^(77)] , 
j 
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where rfi +1 denotes the configuration obtained from r] by interchanging the particles and/or 
vacancies at sites j and j + 1. Following [6U\ I3"T] , attaching a vector space C 3 at each discrete 
point j and using a vector (1,0, 0) for a particle of type A, a vector (0, 1, 0) for a particle of 
type B and a vector (0, 0, 1) T for a vacancy, the operator of the master equation can be written 
as a quantum spin chain which is given by the following expression if one assumes periodic 
boundary conditions 



h =dy: 

3=1 



a<(3 a>/3 a=l 



In the following, the standard vector space 1R will just be called V. In this expression, the 
matrices E" are 3x3 matrices with only one non-zero entry: (E® )^s = dajdps and, as usual, 

L-l 
the expression ^ E c - E-^ means the L-fold tensor product 



ii <g> i 2 <s> ■ ■ • <g> ij-i ® ^ Q/3 ® e^ ® a j+ 2 ® • • • ® 1l- 

The parameters D and q are real and depend on the diffusion rates: q = -v / p^ and D = tJTrTl. 
This Hamiltonian is integrable and the eigenvalues can be found by applying the Bethe ansatz. 
Since the Hamiltonian is non-Hermitian in general, we will encounter complex eigenvalues. 

2.4 Dynamical critical exponent 

In non-equilibrium dynamics, the dynamical critical exponent describes a relation between the 
relaxation time towards equilibrium r (or temporal correlation length) of a system and the 
spatial correlation length £, namely that r ~ £ z with the dynamical critical exponent z. It can 
be shown that for one-dimensional quantum spin chains, r ~ L z . Since the relaxation time is 
dominated by the eigenvalue of the Hamiltonian with the smallest real part (the energy of the 
ground state being equal to zero), we can obtain the exponent z from a finite size analysis of 
the Hamiltonian of the system as 

Re (Ei) = const—. (2) 

L z 

For the single species asymmetric diffusion model, z = I was obtained in [11] from a Bethe 



ansatz calculation for the totally asymmetric case which therefore belongs to the KPZ 
universality class, whereas z = 2 for the partially asymmetric case which describes the Edwards- 
Wilkinson universality class [33]. We will determine the exponent z from a careful study of ([2]). 
We will find the lowest lying eigenvalue of the totally asymmetric diffusion model by means of 
the Bethe ansatz. 

2.5 Nested Bethe ansatz 

We start with the totally asymmetric diffusion model, setting the parameters T^ = and 
Tr = 1. This does not lead to any singularities in the Hamiltonian since its expression contains 
the products Dq = Fr = 1 and Dq~ l = Tl = 0. The new Hamiltonian reads 



" = E 



1 
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-^E(^r^i) 
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This Hamiltonian is integrable, and we will use the algebraic Bethe ansatz to find its spec- 
trum. Although the Bethe ansatz is well known for the Hamiltonian ([I]), the totally asymmetric 
case (given by the Hamiltonian ([3])) cannot be obtained as a special case since setting T^ = 
causes the Bethe ansatz equations to become singular. Therefore, we will derive the Bethe ansatz 
equations for the totally asymmetric case following |34tl35|. We use the following i?-matrix ele- 
ments 

R^a = exp 9 for a = 1, . . . , 3, 

^Ba = ^ s i nn for a < /3, a, /3 = 1 . . . , 3, 

Rf a = fora>/3, a,/3 = 1, . . . , 3, 

R"p = exp 9 for a < j3, a, f3 = 1, . . . , 3, 

R^fa = exp(— 6) for a > /3, a, j3 = 1, . . . , 3, 

where the indices denote the following tensor product in End(V V): 

Ru E m i <g> Ekl, 

or in the language of the associated vertex model this corresponds to the initial state denoted 
by the two indices mi scattering into the final state denoted by the indices kl. 
This i?-matrix satisfies the factorization equation 

R$(e 2 - e 3 )Rf r (e 1 - e 3 )R2 n (ei - e 2 ) = R%{9i - o 2 )s*{o x - e 3 )R%(e 2 - e 3 ). 

We now define the 3 x 3 matrices T(9) , as 

3 

where 

[tab(0)]ij = R%(0), a,b,i,3 = 1,2,3, 
are 3x3 matrices. We define the monodromy matrix as 

/ T [L] T [L] T [L] V A fli ^ 

T [L] = T [L] T [L] T [L] I I Ci Du ^ (4) 

V T§ T§ T§ ) \C 2 D 21 D 22 J 

The monodromy matrix satisfies the fundamental relation 

R{6 - 9')[T^(9) ®tM(9')] = [TW(0') ®tW{0)]R(8 - 9'). (5) 

The transfer matrix 
3 

r(9) = Y J Td9) 

i=l 

is obtained as the trace of the monodromy matrix. It is a matrix acting on V 8 . In the above 
notation, it can be written as 

t(8) = A(8) + D 11 (9) + D 22 (8). 

We recover the Hamiltonian ([3]) as the logarithmic derivative of the transfer matrix, if the 
spectral parameter 9 is set to zero: 

rf(logT(fl)) 

d9 
Therefore, diagonalizing H will be the same problem as diagonalizing r. 



H 



(6) 
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2.6 Diagonalization of the transfer matrix 

We start be deriving algebraic relations for the matrices A, Bj, Cj, j = 1, 2 and Dik, i, k = 1, 2 
in equation (j4|). 

Since we would like to diagonalize r, we are looking for a reference state that is a simultaneous 
eigenstate of A and Da and is annihilated by C% and Dij for i ^ j. 

We choose 





O) = ®0 ® • • • ® = ®fc=i (7) 



to be our reference state. This is not necessarily the physical ground state. T^ L ' acting on this 
reference state is interpreted as follows: 

Tg ] \n) r^|fi) T^\n) \ / a (9) L \n) B^n) B 2 \n) \ 
T^\n)= rg ] |fi) r|J|o) rg ] |o) = o c(e) L \n) o , (8) 
rgi|n> rj?|n> rj?|n> / V ° ° C ( 6, ) L I«) / 



where a(0) = exp(#) and c(8) = 2sinh(0). The operators Bi(9,a) will create new states when 
acting on the reference state. They will be called creation operators in the following. 

In the next step we use the fundamental relation equation ([5]) to derive the bilinear algebra 
among the operators A(8), Bi(8), Ci(9) and Dij(8). 

The relations we are going to use later on are: 

A(9)B i (8') = g(8' - 8)B i {8')A{8) - h{8' - 8)B i (8)A(8 f ), 

B i (p)B j {0) = r%B p {0')B q (e), 

D i:j {9)B k {d') = g{8 - 8')rf q B p (8')D n (8) - g{8 - 0')B i (e)D jk {0'). (9) 

Here r % p q are coefficients of a 4 x 4 matrix where only the following five entries are different from 
zero: r\\ = 1, rjf = 1, rfl = 2smh(8)e~ e , r\\ = e~ 2e , r|| = 1 and the functions g(8) and h(8) 
are given by: 

e e i 

g (8) = , h{8) = 2e- e sinh(0) 



2sinh(<9)' w w g{8)' 

Now we are ready to look for an eigenstate of the transfer matrix t{8). 

We will make the following ansatz for an eigenfunction of the transfer matrix having the 
creation operators Bi(8), i = 1, 2 act on the reference state O given by equation ([7]). 

^(Ai,A 2 ,... , Xp) = /^ XffB ff ( yi ) (Ai)-Bo.fa 2 ) (A2) ■ ■ ■ B a (y p )(\ p )\n) , 
M 

where the sum runs over all possible permutations of the set of indices {2/1,2/2, • • • , Up} £ {1)2}. 
We will calculate the action of the transfer matrix r on the vector $. Recall that r = 
A(8) + Dn(8) + -£>22(#)- So we will encounter terms of the form 

A(8)B a{yi) (X 1 )B a{y2) (X 2 ) ■ ■ ■ B a(yp) (X p )\n) and 
Ai(^)-B ( r( 2 /i)(Ai)-B -(j /2 )(A2) • • • B a (y p )(X p )\n). 

We know how A{9) and Da(8) commute with the operators Bi(Xj) from the relations ([9]) and 
we know how A(8) and Da{8) act on £1 from equation (J5J). 
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We will get two types of terms. In the first type of terms, we get the original combination 
of B a /y 1 \{Xi)B a r y2 ^{X2) ■ ■ ■ B a ( y )(A P )|0) back. They originate from taking the first term in the 
commutation relation for Bi(9)Bj(9) in Q. Using standard terminology, we will call those 
terms wanted terms. In the second type of terms, one of the B a i yi \ operators will depend on the 
parameter 8. These terms will be called unwanted terms and their sum has to be set to zero. 

We will also use the following notation: We define a 2 P dimensional vector B containing all 
possible terms of the form B a ^ yi ^{X{)B (T ^ y2 ^(X 2 ) ■ ■ ■ B ff ( y ) where, as above, a denotes a permu- 
tation of the set of indices {y±,y 2 , ■ ■ ■ , y P } £ {1,2}. We also define a 2 P dimensional vector X 
containing the corresponding x a in the same order. This allows us to rewrite the sum as a scalar 
product of these two vectors: 

Y,^B a(yi) {\ 1 )B a{y2) (\ 2 ) ■ ■ ■ B a{yp) {X p ) = £(Ai) ® B{\ 2 ) ® ■ ■ ■ <g> B{X P ) -X = M-X. 
M 

Then the action of the transfer matrix t{9) on ^(Ai, A2, ■ ■ • , X p ) leads to wanted terms that 
can be written as: 

v 

a(e) L l[g(x l -e)^(x 1 ,x 2 ,...,x p ). 

1=1 

The unwanted terms arise from taking the second term in the commutation relation for the 
B{9) operators ([9]). E.g. one of the unwanted terms appears if the second term in the commutator 
is applied to A(6)B(X k ) and is of the form 

v 
h(X k -6)a(X k ) L ~[[g(X n -X k ). 

As observed by de Vega [36] the general unwanted term can easily written down if one uses 
the following fact about a cyclic permutation of p operators B(Xj): 

fl(Ai) ® B{X 2 ) ® • • • ® B(X P ) = B(X 2 ) ® B(X 3 ) ® • • • ® B{X P ) ® B^XMJXi, {A}), 



2 



M, 



where rf 2) (#,{A}) = £ 7^(0, {A}) and 

a=l 

r&ty, {A}) = £ 41 (* - Ai) ® 41,^ - A 2 ) ® ■ ■ ■ ® 4t l6 (0 - A p ), 

ai,...,ap_i=l 

[21 

where [iL(^)]y, a ,b,i,j = 1,2 are 2x2 matrices where only the following five entries are 
non-zero: 

[t?l(e)]u = rll(9), [tf}(6)]22 = rg(0), [tjg(0)]ai = rg(0), 
[*SW]ia = rgW, [t|W]22 = rI(0). 

Here r^AO) are the coefficients appearing in equation and t,^ (#, {A}) is the transfer matrix 
of the six- vertex model (i.e. a model with two states instead of three states) for a line of p sites 
with inhomogeneities {Aj}, i = 1, . . . ,p. 

So the cyclic permutation B(Xi) — > B(Xi+i) followed by the multiplication of the matrix 

M = r ( p 2) (A 1 ,{A}) 
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leaves the product B(X±) ® -B(A 2 ) ® • • ■ ® B(X P ) invariant. Therefore the general unwanted term 
in the result for t4(0) x I'(Ai • • • X p ) can be written as 

p v 

-Y,{K^k-e)a{\ k ) L J] g{\ n - X k )B{9) ® J3(A fc+ i) ® ■ ■ ■ 

A;=l n=l,n^k 

® B(Ap) £(Ai) ® • • • ® ^(Afc.ijM^ 1 }!^}. 

Similarly one can find the wanted and unwanted terms after acting with (Dn(9) + 1)22(0)) 
on |\&). The wanted term is 

p 

(D n (e) + r>22W)l*> = Ug(o - a 3 -)c(0) l b(Ai) ® b(a 2 ) ® • • • ® b(\ p )t$(9, {\})x\n). 

Using the same argument as before, the general unwanted term can be written as a sum of terms 
where B{9) replaces B(X k ): 

p p 

£{-M0-Afc) II g{^k-\n)c{\ k ) L B{6)®B{\ k+l )®--- 

k=l n—=l,nj^k 

® B(Ap) £(Ai) ® • • • ® ^(Afc-OM^V^CAfc, {A})}A|ft>. 

Altogether, the sum of wanted terms reads 

p p 

S(Ai) B(A 2 ) • • • ® B(A p )[a(0) L [J 5 (Ai - 9) + \[g{9 - X j )c(9) L r^(9, {\})]X\Q). 

i=i 3=1 

The sum of unwanted terms reads: 
p 



^2[B(0) <g> B(X k+l ) ® • • • ® B(X P ) ® B(Ai) ® ■ ■ ■ ® B(A fe _!)]M* 



x J /i(A fe - #)a(A fc ) L J] 9(K - X k ) + h{9 - X k ) J] 5 (A fc - A n ) 

V n=l.nj^k n—=l,nj^k 

xc(x k ) L r$(x k ,{X})\x\n). 

Since we are looking for an eigenstate of the transfer matrix, the sum of all wanted terms has 
to be proportional to \^f) and the sum of all unwanted terms has to be equal to zero. The first 
condition determines X to be an eigenvector of r^ corresponding to the eigenvalue A( 2 )(#, {A}): 

t$X = A {2) (9,{X})X. (10) 

The eigenvalue A/ 2 \(9, {A}) is determined by requiring the sum of unwanted terms to become 
zero: 

We have reduced the original eigenvalue problem to equation §W\\ which is an eigenvalue equation 
for the transfer matrix tvIV of the six-vertex model. This is the crucial step in this paper. All 
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we have to do now is repeat the same diagonalization procedure for equation (|10p . We repeat 
exactly the same steps as before, acting with the transfer matrix on a reference state with r 
creation matrices B, and deriving equations for the wanted and unwanted terms. 

This leads to another expression for the eigenvalue A( 2 ) that can be set equal to equation (jlip 
and another consistency equation stemming from setting the unwanted terms equal to zero. 
In this way we obtain the following nested Bethe ansatz equations which have two types of 
unknowns, A&, k = 1, . . . ,p and Aj, j = 1, . . . , r, 



exp(Afc) 



2sinh(A/ ) 



n<-«p<«.-«.»n^& 



k=l,...,p, 



fc=l J n—l 

The eigenvalues of H are obtained using the relation between the transfer matrix and the 
Hamiltonian given by equation (jBJ): 

* = SMJWWU = * + E3ffii&j- < 13 > 

2.7 Numerical solution of the Bethe ansatz equations 

Since the Hamiltonian ([3]) is describing a reaction-diffusion system, the ground state is the 
steady state with energy zero and all excitation energies will have positive real parts. We 
fixed L, p and r and solved the coupled system of Bethe ansatz equations numerically. We also 
diagonalized the Hamiltonian ([3|) numerically for a small number of sites (typically up to L = 9) . 
We compared the energies obtained from the Bethe ansatz with the ones obtained by numerical 
diagonalization of Hamiltonian H to single out the eigenvalue with the smallest non-zero real 
part which plays the role of the second lowest eigenvalue and determines the gap. This first 
excited state always lies in the sector with p = ~, r = 0. This state has an equal density of 
particles A, B and vacancies of ^. We extrapolated the energy values of the first excited state 
for L —?• oo. 

2.8 Results 

We found the solution for p = -g, r = corresponding to the first excited state of Hamiltonian. 
To make sure that this really is the lowest lying excited state, we compared to the results of 
a numerical diagonalization of the Hamiltonian for finite lattice lengths. This state has an equal 
density of particles A, B and vacancies of 3. The roots all lie in the complex plane. A typical 
pattern of roots for the first excited state of the Hamiltonian is shown for L = 36 in Fig. [2j The 
solid line is the unit circle that is given as a guide to the eye. 

The data was extrapolated using the Bulirsch-Stoer algorithm [37] . Since we expect a scaling 
of the form (remember the ground state energy is always zero) 

Re(AE(L)) = Re{E 1 (L)) = const L~ z + o(L~ z ), 

we built extrapolants for the exponent — z as 

h °&[Re(AE(L+3)) 
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Figure 2. Complex roots Afc, k = 1, . . . , 12 for the Bethe ansatz equations with p = j, r = and 
L = 36. 

The data used in the extrapolation for the exponent is shown in table: 

L extrapolant 

6 -1.6336892192762 

9 -1.6252314332778 

12 -1.6092183117219 

15 -1.5952666540982 

18 -1.5839870664789 

21 -1.5749003909369 

24 -1.5674968193872 

27 -1.5613778750522 

30 -1.5562495252464 

33 -1.5518961566109 



The result of the extrapolation with BST-algorithm is -z = -1.50000009 with error 0.00000323. 
This clearly shows that the exponent z is |. 

It is interesting to note that this state corresponds to the Bethe ansatz equations with only 
one type of roots. The vanishing of the second type of roots does not correspond to the simple 
inclusion of the one-species model into the two-species model given by making the density of the 
second type of particles zero but rather for the state we consider corresponds to equal densities 
of these in-equivalent particles. The fact that the densities are coupled may suggest that there 
is some underlying quasi-particle formalism which may give a theoretical way to explain the 
connection to the single species exclusion model and the occurrence of the exponent | well 
known from the totally asymmetric single species exclusion model. 



2.9 Analytical solution of the Bethe ansatz equations 

First we would like to rewrite the Bethe ansatz equations equations (|12p and (|13p in integral 
form. To that end we use the following two changes of variables: e k = Zk, z\ = Z&, k = l,...,p 
and e A J = y,-, y? = Yj, j = 1, . . . , r. Equation (fT2|) becomes: 



Vj, Vj 



Zi- 



Z k -l 



m-tin 



3 = 1 



Y, 



3=1 



(Yj - Z k 



fc = l. 



,P, 
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Figure 3. Transformed complex roots Zk, k = 1, . . . , 12 for the Bcthe ansatz equations with p = jf, 
r = and L = 36. 



n 



Y 



Y - Z k 



m-% 



n+3 



(14) 



The new equation for the energies reads: 



E = L + ^ 



fc=i 



2Z k 
Z fc -1' 



Our nunrerical work described above suggests that the first excited state lies in the sector with 



3 ' 



0. After applying the above transformations, the transformed roots Z k , k = 1, . 



L 
' 3 



V 

will lie on a curve that is shown for L = 36 in Fig. [3l 

In order to analyze the Bethe ansatz equations in the limit of large lattice lengths L, it is 
convenient to introduce a function 



g(z) = In 
and a function 
K(zi,z) = 



z-\ 



lnl- 



2/ 



For both definitions the branch cut of the ln-function is taken along the negative real axis. The 



Bethe ansatz equations for p 



3 ' 



can now be written as 



YUZj) 



L J ' 



3 = 1, • • ■ 



L 



with a so-called counting function 

L/3 

iy L (z) = 5 (z) + -^K(^,, 



/=i 



Each excited state of the Hamiltonian in the sector with p 
set of integers {Ij \ j = 1, . . . , g- }. 



will correspond to one 
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Figure 4. Sketch of the integration contour. The open dots represent the roots Zj 



The situation is very similar to the one described in [38] for the one-species asymmetric 
exclusion model. Therefore we can adopt the technique to transform the discrete Bethe ansatz 
equations (|14p to integral equations by using an identity that follows from the residue theorem: 



i L/3 i r /1 

rEM) = J" f dzf(z)Y[(z)cot(-LY L ( i 



In this equation, C is the contour enclosing all the roots Zj. We can view this contour C as 
the union of two contours C\ and Ci as shown in Fig. [H 

These two contours intersect at two points £ and £*. We will fix those two points by requiring 

Y L {e) = -*+\, yuo=k-1- 

Rewriting this integral by separating the contributions coming from two contours C\ and C2, 
this becomes: 



1 l K(v,,,)YL 



1 ft If K(w z)Y' 

iY L {z) = g(z) + -^K( W ,z)Y[( w ) dw + -^ -Ij^- cto + ^^(^ 

The formula for the energy reads: 



~ / e(z)Y[(w)dw + " 



2tt 



£* 



<*)?L 



dw + 



L 



2vr / ri 1 - e -iLY L (w) 2vr Jr„ e iLY ^ w ) - 1 



Ci 



*(*)*z 



-diu, 



C 2 



where the function e(z) is given by 

2z 
e(z) 



z-1 



It should be possible to analyze these equations along the lines of de Gier and Essler |38j by 
an expansion in inverse powers of L; this will be left for a future publication. 
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3 Conclusion and outlook 

We have shown that the dynamical critical exponent of the totally asymmetric exclusion model 
is 2 which in the single-species asymmetric diffusion model is the exponent for the KPZ uni- 
versality class. This is a new result that cannot be deducted from the recent paper by Arita et 
al. [25]. 

The Bethe ansatz equations and their solutions are qualitatively very different from the 
case of a single species asymmetric exclusion model. It would be very interesting to derive 
this result analytically once the patterns of the solutions of the Bethe ansatz equations are 
fully understood. Building on the results reported here, the next step would be to solve the 
integral equations representing the Bethe ansatz equations in the continuum along the lines of 
de Gier and Essler [38]. Since the integration is performed over the curve formed by the roots 
in the thermodynamic limit, a qualitative understanding of this curve is a prerequisite for this 
calculation. We hope to report on this soon. 

Another new direction will be to generalize the Bethe ansatz to the case of open boundaries. 
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